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Abstract. As the first gravity wave (GW) climatology study 
using nadir-viewing infrared sounders, 50 Atmospheric In- 
frared Sounder (AIRS) radiance channels are selected to es- 
timate GW variances at pressure levels between 2-100 hPa. 
The GW variance for each scan in the cross-track direction 
is derived from radiance perturbations in the scan, indepen- 
dently of adjacent scans along the orbit. Since the scanning 
swaths are perpendicular to the satellite orbits, which are in- 
clined meridionally at most latitudes, the zonal component 
of GW propagation can be inferred by differencing the vari- 
ances derived between the westmost and the eastmost view- 
ing angles. 

Consistent with previous GW studies using various satel- 
lite instruments, monthly mean AIRS variance shows large 
enhancements over meridionally oriented mountain ranges 
as well as some islands at winter hemisphere high latitudes. 
Enhanced wave activities are also found above tropical deep 
convective regions. GWs prefer to propagate westward above 
mountain ranges, and eastward above deep convection. AIRS 
90 field-of-views (FOVs), ranging from +48° to —48° off 
nadir, can detect large-amplitude GWs with a phase velocity 
propagating preferentially at steep angles (e.g., those from 
orographic and convective sources). The annual cycle domi- 
nates the GW variances and the preferred propagation direc- 
tions for all latitudes. Indication of a weak two-year varia- 
tion in the tropics is found, which is presumably related to 
the Quasi-biennial oscillation (QBO). 

AIRS geometry makes its out-tracks capable of detecting 
GWs with vertical wavelengths substantially shorter than the 
thickness of instrument weighting functions. The novel dis- 
covery of AIRS capability of observing shallow inertia GWs 


will expand the potential of satellite GW remote sensing and 
provide further constraints on the GW drag parameterization 
schemes in the general circulation models (GCMs). 


1 Introduction 

Gravity waves (GWs) are known to play a key role in global 
climate and weather dynamics by transporting energy and 
momentum from the lower to the upper atmosphere, trans- 
port which is essential for determining the general circula- 
tion and temperature structure in the stratosphere and meso- 
sphere (e.g., mesosphere wind reversal, quasi-biennial oscil- 
lation, etc.). GWs are of particular importance for wave dy- 
namics in the summer hemisphere in the stratosphere where 
planetary waves are weak, and are important in the entire 
mesosphere. As an example, model simulation results sug- 
gest that the mesosphere summer easterlies would arrive one 
month late without including gravity wave drag (GWD) in 
the model (Scaife et ah, 2002). GWs also affect weather and 
chemistry. For example, temperature fluctuations associated 
with GWs can lead to formation of polar stratospheric clouds 
(PSCs) in a synoptically warm condition and susequently af- 
fect ozone depletion (Hamill and Toon, 1991). 

The importance of gravity waves on climate and weather 
became widely recognized in recent decades due to progress 
in observational techniques and model improvements. Two 
issues remain as major concerns, however. Firstly, a signif- 
icant amount of the GWs cannot be resolved in GCMs and 
hence need to be parameterized except those with long wave- 
lengths and low frequencies. These prescribed parameters, 
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formulated to be dependent on detailed characteristics of 
GWs and GW sources, however, are poorly constrained by 
observations, and are actually heavily “tuned” in most of the 
cases to attain a realistic atmosphere. For example, one of the 
major GW sources, convection, is often assumed to be uni- 
formly distributed all over the globe in GCMs with spectral 
GWD parameterization schemes, which is obviously far from 
reality (Alexander and Rosenlof, 2003). The other important 
non-stationary GW source - jet imbalances - is not included 
in most of the GCMs (Kim et al., 2003). Secondly, different 
observational instruments can only see partial gravity wave 
spectra. Although solving this jigsaw puzzle is underway 
(e.g., Alexander et al., 2010), the high-frequency portion of 
the spectrum is invisible to most of the instruments (Wu et al., 
2006; Alexander et al., 2010), and the momentum flux asso- 
ciated with these waves remain uncertain. GWs of higher in- 
trinsic frequency and shorter horizontal wavelengths can po- 
tentially carry greater momentum flux than lower frequency 
GWs, and could therefore exert significant wave drag (Fritts 
and Alexander, 2003). Currently we are missing global infor- 
mation in particular on these shorter horizontal scale waves. 
This presents an important information gap for any effort to 
constrain gravity wave drag (GWD) in climate models by 
measurements. 

Among various GW observations, satellite instruments 
have advantages over others (e.g., radiosonde, radar, aircraft 
measurements) in terms of high and regular sampling den- 
sity and global coverage. As an example, the Aqua satel- 
lite, which carries the Atmospheric Infrared Sounder (AIRS), 
samples each 2° x 2° latitude-longitude gridbox at least 24 
times per month. Intercomparisons and cross-validations 
can therefore be readily carried out among different satel- 
lite instruments on a global basis. Based on their view- 
ing geometry, satellite GW observational instruments can be 
further categorized into three groups: limb, sub-limb, and 
nadir sounders. The limb sounders (e.g., LIMS, CRISTA, 
HIRDLS, GPS) are sensitive to GWs with A. z /A/j ratio of 
less than 1/10-1/20, i.e., low-frequency GWs. Here X z and 
X x denote vertical and horizontal wavelength, respectively. 
The nadir sounders (e.g., AIRS, AMSU-A, SSMIS) are most 
sensitive to high-frequency GWs, and the sub-limb sounders 
(e.g., MLS) are sensitive to mid-frequency GWs. Table Al 
summarizes the sensitivities of wave parameters from four 
satellite instruments that are compared with AIRS results 
later on. More details can be found in Wu et al. (2006) and 
Preusse et al. (2009). 

As a nadir viewing sounder, AIRS has at least three supurb 
properties for GW studies. Firstly, AIRS has a high horizon- 
tal resolution (~13 km at nadir). Therefore AIRS is very sen- 
sitive to high-frequency GWs, which are badly constrained in 
their global distributions by both observations and models. 
Secondly, for the outer tracks the sensitivity depends on the 
tilt of the wave fronts. This allows us to estimate preferen- 
tial GW propagation directions from the viewing-dependent 
variance difference between the two outmost off-nadir views. 


Lastly, AIRS radiances are a direct measure of GW induced 
air temperature perturbations, and hence provide a more ac- 
curate measurement of GWs than retrieved temperatures. 

This paper is organized as follows. The following section 
gives a brief introduction of the AIRS instrument, and a de- 
scription of the methods we use to retrieve GW properties. 
In Sect. 3, a climatology of wave variance and the preferred 
zonal propagation direction will be derived in terms of zonal- 
mean, geographical distribution and temporal variations. In- 
terpretation of the results is discussed in detail in Sect. 4, 
where we emphasize on two major gravity wave sources - 
topography and convection, and their interannual variabili- 
ties. Wave enhancements in the equatorial lower stratosphere 
and the winter pole at lOhPa will be discussed in a great 
length, and a numerical experiment is performed to support 
the speculation on the observed GWs being low-frequency 
waves. Concluding remarks can be found in Sect. 5. 

For reader’s convenience, we hereafter refer to “inter- 
nal GWs” as high-frequency (a>^>/), long vertical wave- 
length (/,,> 1 0 km) and short horizontal wavelength waves 
(A/,~100km or less), and use “inertia GWs” to rep- 
resent low-frequency (a>~/), short vertical wavelength 
(a- < 1 0 km) and long horizontal wavelength (A/, ~ 1000 km) 
waves, where co and / denote wave frequency and the Cori- 
olis parameter. One should keep in mind that GWs in the 
stratified atmosphere are essentially all internal waves. 

2 AIRS instrument and method 

AIRS is an infrared spectrometer and sounder that con- 
tains 2378 channels in 3.74-4.61 pm, 6.20-8.22 pm and 8.8- 
15.4 pm wavebands. It makes a cross-track scan every 8/3 
second that includes 90 footprints on the ground, 4 indepen- 
dent views over the cold space, and 3 views into 3 differ- 
ent calibrators (Aumann and Miller, 1994). The scan swath 
is ~1600km wide, reaching ±48.95° from nadir. The an- 
gle difference between two adjacent footprints, therefore, is 
about 1.1°, which corresponds to ~ 13 km footprint size at 
nadir. The satellite orbits over 24 h are divided into 240 gran- 
ules, each of which contains 135 scans (6 min duration). 

We use AIRS CO 2 15 pm radiance emission bands with 
wavenumbers ranging from 665 cm - 1 to 693 cm - 1 in this 
study. 50 channels are selected with 1 1 distinguished weight- 
ing functions (WFs) peaking between 2hPa and 100 hPa. 
The channel numbers are listed in Table A2 in the Ap- 
pendix A. Ideally, only channels at the maxima or minima 
of the radiance spectra, or paired channels at two side-wings 
of extremes should be used since they are most stable with 
small detector drifts on the spectra, which mainly occur 
due to the expansion/contraction of the glass filter depend- 
ing on whether it faces the Sun or in darkness (i.e., ascend- 
ing/descending of the orbits). However, we are losing infor- 
mation over most of the middle to upper stratosphere by only 
selecting these channels. With this consideration, several 
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channels on one side wing of maxima/minima (bold num- 
bers in Table A2) are also selected. Results from ascending 
and descending orbits have been compared carefully, and no 
significant differences have been found whatsoever, so this 
problem shouldn’t cause major problems or failures of using 
the wing channels in our research. 

Since AIRS scans are always perpendicular to the orbit 
track, the scan line can therefore be treated approximately 
as along the west-east direction for most of the time, except 
at high latitudes where the scan is mostly meridional (Aqua 
satellite orbit track can be checked at http://www.ssec.wisc. 
edu/datacenter/aqua/). We hence exclude the data beyond 
±80° latitudes. 

The variance of the radiance (cr 2 ) has three components: 

^^GW+^noise + f 2 ( ! ) 

Normally other variance (c 2 ) is much smaller than the first 
terms on the right-hand of Eq. (1), and hence is omitted here. 
Evidence shows that AIRS may observe some turbulence sig- 
nals or small vertical wavelength inertia GWs, which will be 
discussed later on in Sect. 4.3. However, these additional 
sources of variance are highly uncertain, and we would rather 
include them in rr ( 2 w than separating them out. Since the in- 
strument noise is on the same order with the magnitude of 
the GWs, it needs to be carefully evaluated and subtracted 
consequently. The instrument noise (cr 2 oise ) is estimated in 
an analagous way to Wu and Waters (1996). Instead of us- 
ing pre-launch information, we feel it is more accurate to 
estimate the random component of the noise from the real 
data with the assumption that cr 2 oise only depends on the fre- 
quency channel, and should not vary with location, time and 
viewing angle. To extract the radiance noise, we first apply 
a 3-pt running smooth window to each independent scan and 
compute the variance as the squared difference between the 
original and smoothed data (referred as “3-pt result” here- 
after). Then, the minimum variance from the monthly aver- 
aged 2° x 2° map for each month is determined. Because the 
3-pt variance is affected least by GWs, the minimum vari- 
ance is assumed to be very close to the measurement noise. 
These minimum variances are found always coming from 
the summer hemisphere high-latitudes. Finally, the mean is 
computed for all the variances that are below the median of 
the timeseries, and we attribute this mean as the instrument 
noise cr 2 oise . The 3-pt variance is used for the noise estima- 
tion since it gives the smallest cut-off wavelength and hence 
contains the most high-frequency information. Nevertheless, 
the 3-pt result still contains a fair amount of GW contribu- 
tions in that the latitudinal and geographical distributions of 
the 3-pt result look similar to 7-pt result even without noise 
subtraction (not shown). The noise estimation method used 
here is confirmed to be reasonable and accurate by the fact 
that variance reaches minimum at summer hemisphere high 
latitudes, and is quite uniformly and stably distributed there 
throughout the time. Moreover, there are no statistically sig- 
nificant differences among different viewing angles there. 


We furthermore compared the noises estimated by the 
above method with the NEdT (Noise Equivalent Delta Tem- 
perature) listed in AIRS property reports (http://airs.jpl.nasa. 
gov/data_products/algorithms/), and as one can easily tell 
from Table A2 in the Appendix A, they agree very well. 
A linear increasing trend has been found in the estimated 
noise at various levels, which is on the same order of the 
trend found in the AIRS documentation. Since this trend is 
far from significant (regression coefficients on the order of 
10 -5 K 2 month -1 ), which is probably due to the run-off of 
the sensors, we do not remove the trend in this study. 

As suggested by Wu and Waters (1997) and Wu and Eck- 
ermann (2008), the minimum detectable GW variance (a ( 2 w ) 
is a^ m — «J2 (M — 2) //Vcr 2 oise , where M is the number of 
points used in this filter for variance estimation (M — 1 in 
this case). Since cr 2 oist , is determined from a monthly mean 
2° x 2° map, a zonal average with N points in this grid size 
can further reduce the detection threshold. The outmost M/2 
points are used to get an averaged variance value, which 
also needs to be considered. The minimum detectable val- 
ues at various pressure levels are listed in Table A2 in the 
Appendix A. 

We cannot directly calculate the GW horizontal propaga- 
tion direction from radiance measurements of a single scan. 
However, the preferred zonal propagation direction can be 
indirectly inferred from the difference of GW variances be- 
tween two different FOVs. This is illustrated in a schematic 
picture Fig. 1, where a typical mountain gravity wave (MW) 
packet is generated downstream of the mountain under a 
westerly wind. Wave phases, though stationary in the Earth 
frame, are propagating downward and westward relative to 
the westerly mean flow, while energy is propagating upward 
and westward in this case. In the west field-of-view (FOV) 
of AIRS, the wave amplitudes are largely smeared out, while 
the east FOV would yield a much greater reading. There- 
fore, negative (positive) value of (GW w — GW e ) reflects a 
predominant westward (eastward) zonal propagation. Since 
multiple factors affect the magnitude of the difference, such 
as the viewing angle relative to the phase front, the original 
GW amplitude, etc., the absolute value of the difference does 
not have a definitive physical meaning. As long as it is sta- 
tistically significant, the sign of the difference is the variable 
that really matters. This approach has been formerly used on 
Aura MLS observations to study GW meridional propaga- 
tion direction, and has turned to be a very useful and effective 
method (Wu and Eckermann, 2008). 

Before presenting the results, one should note that the 
majority of the GW signals in AIRS comes from high- 
frequency GWs. The visibility can be computed as a func- 
tion of along-track wavelength (X x ) and vertical wavelength 
(A-) by convolving a plane wave of unity amplitude (IK) 
with the weighting function (WF) (McLandress et al., 2000). 
In the along-track direction, the peak of the sensitivity is 
mainly determined by the length of the truncation window, as 
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West FOV East FOV 

Fig. 1. Cartoons showing the wave smearing effect from different 
views, (a) Westmost off-nadir view, where in this particular case 
wave signals are largely smeared out; (b) Eastmost off-nadir view, 
where the strongest variance should be observed. Red parallel rect- 
angles represent the phase lines of the GW, while bold red arrows 
denote the mean wind direction. 


suggested by Fig. 2. The left half of the three curves, which 
is similar for the 3-pt, 7-pt and 15-pt results, is controlled 
by the width of the AIRS FOV. All three rise sharply once 
X x becomes larger than the footprint size, which is ~ 13 km 
at nadir. The length of the smoothing window decides the 
peak as well as the tail of the filter. In the vertical direc- 
tion, AIRS can hardly detect GWs with X z shorter than the 
full width at half maximum (FWHM) of the WFs, which is 
~ 12 km. However, there are some exceptions, which will be 
discussed later in Sect. 4.3. The sensitivity increases mono- 
tonically with the increase of the X z before reaching the peak 
at X z ~ X x max , where A* max denotes the most sensitive along- 
track wavelength. This can be better understood from revis- 
iting Fig. 1, where the largest amplitude would be reached 
when the GW fronts are parallel to the viewing angle. For a 
7-pt window, the peak sensitivity occurs when both X x and 
X z equal to ~ 100 km. Alexander and Barnet (2007) stud- 
ied multiple MW events over the tip of Andes and Antarctic 
Peninsula, and found out the most dominant X x is around 
100-150 km. So the 7-pt window can effectively capture 
these GWs. With a wider window, we are at the cost of pos- 
sible inclusion of other waves (e.g., Kelvin wave), and we 
may lose the dominant GWs as suggested by Alexander and 
Barnet (2007). Nevertheless, the major conclusion’s robust- 
ness does not depend much on window size. A 3-pt analyses 
figure is provided in the Appendix B, which shows smaller 
amplitude and similar patterns to the 7-pt result, implying 
that 7-pt is a better window to use. We hence only present 
7-pt results in this paper. 


3 Climatology of gravity wave variances and the 
preferred zonal propagation direction 

The monthly mean latitudinal, geographical and temporal 
variations will be presented in this section. Year 2005 is 
used to indicate the mean GW climatology, which represents 



HORIZONTAL WAVELENGTH (KM) 


Fig. 2. AIRS visibility as a function of the along-track wavelength 
with 3-pt (black), 7-pt (blue) and 15-pt (red) running smooth win- 
dow applied. 


a condition with no strong sudden stratospheric warming nor 
ENSO event. 

3.1 Latitudinal distributions 

Figures 3 and 4 plot the monthly averaged zonal mean 7- 
pt radiance variance as a function of latitude and height in 
January and July of 2005, respectively (with estimated in- 
strumental noise deducted), overplotted with UK Met Of- 
fice monthly mean zonal winds. The variances can be ex- 
plained by GWs, since they in general grow exponentially 
with height at mid-high latitudes. In the middle to upper 
stratosphere (above 20 hPa), the two largest variance centers 
colocate with the westerly and easterly jet cores in the up- 
per stratosphere, and the one in the winter hemisphere is ap- 
parently stronger. This indicates that large amplitude GWs 
(or frequently generated GWs) are from sources beneath the 
two jet centers. The overall feature is similar to previous 
GW studies using Aqua AMSU-A, UARS MLS and Aura 
MLS (Wu et al., 2006) with comparable strengths in the win- 
ter hemisphere, but weaker amplitudes in the summer hemi- 
sphere. However, the wave enhancement at the lOhPa near 
the winter pole is a new feature revealed by AIRS data. An 
interpretation of this feature is inertia GWs, which will be 
discussed in Sect. 4.3. Very low GW activity is observed by 
AIRS in the summer polar regions. Even without the sub- 
traction of estimated noise, the variance there is still signifi- 
cantly smaller than GW variance derived from other satellite 
measurements (Wu and Eckermann, 2008). Low phase speed 
GWs cannot propagate into these regions because of the 
wind reversal between the troposphere and the stratosphere. 
In addition, the low wind velocities in the summer high- 
latitude stratosphere allow AIRS only to observe GWs of 
high ground-based phase speeds (larger than ~40 ms -1 ; cf. 
discussion of observational filter by e.g. Alexander, 1998). 
However, such fast waves should be emitted from sources 
such as fronts that are currently parameterized in models to 
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Fig. 3. Zonal mean of the GW variance at the westmost view (a), eastmost view (b) and the difference between the westmost and eastmost 
views (e) as a function of latitude and height for January, 2005. Monthly mean zonal winds obtained from UK Met Office are contoured in 
solid (westerly) and dotted (easterly) lines with an interval of 10 m s _ 1 . The color scale is linear by taking a square-root of the values. In (c), 
positive (negative) values correspond to preferred eastward (westward) propagation direction. 


(o) West (b) Eost (c) West - Eosl 
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Fig. 4. Same with Fig. 3, except for July 2005. 


reproduce the wind reversal around the mesopause (Richter 
et ah, 2010). The very low values in AIRS indicate rather low 
probability of such waves to occur and support observational 
evidence that GWs from subtropical convection propagating 
polewards are a major source of momentum for the high lat- 
itude summer mesosphere (Preusse et ah, 2009). 

In the lower stratosphere, AIRS GW variance is generally 
larger in the tropics and decreases with latitude. Similar re- 
sults have been found in some measurements that represent 


the behavior of the inertia GWs (e.g., Alexander et ah, 2002; 
Wang and Geller, 2003; Ratnam et ah, 2004; Preusse et ah, 
2006) and middle frequency GWs (e.g., Wu and Eckermann, 
2008). AMSU-A, which has the same scan angle as AIRS, 
and hence which is supposed to be also sensitive to internal 
GWs, observes features opposite to AIRS at the equatorial 
lower stratosphere where the GW variance reaches minimum 
(Wu et ah, 2006). This interesting feature will be discussed 
in details later on in Sect. 4.3. 
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(d) West, 10 hPo 






(c) West - Cost, 2.5 hPo x io-’(k’) 



(f) West - East, 10 hPa X 10“ (K’) 



(i) West - East, 80 hPa x 10“ (K‘) 




Fig. 5. Geographical distribution of GW variance on a 2° x 2° grid at westmost view (left column), eastmost view (middle column) and 
the difference (right column) for January 2005 at 2.5, 10, and 80hPa from top to bottom. 3-pt smoothing has been further applied to make 
the signals stand out. In the left two columns, values smaller than 2cr are uncolored, where a is the minimum detectable variance. In the 
rightmost column, absolute values smaller than the minimum detectable difference are whitened. The UK Met Office monthly mean zonal 
winds at corresponding levels are contoured in black solid (westward) and dotted (eastward) lines. Total wind velocity greater than lOms^ 1 
areas at 700 hPa are hatched with brown lines. 


As stated in Sect. 2, one can infer the GW preferred 
zonal propagation direction from the difference of the GW 
variances between the two outmost views. As in the right 
columns of Figs. 3 and 4, the zonal mean zonal propaga- 
tion directions, on a monthly scale, always tend to be in the 
opposite direction to the mean zonal winds, which is in gen- 
eral westward (eastward) in the winter (summer) hemisphere. 
This is true especially for large amplitude GWs at high alti- 
tudes, as will be shown in the geographical patterns in the 
next section. Orography is one of the most prominent GW 
sources in the winter hemisphere. Only in strong westerly 
winds can mountain waves, which have zero ground- based 
phase speeds, gain high intrinsic phase-speeds and long ver- 
tical wavelengths, and thus become visible to AIRS. Notice 
that the overall pattern in Fig. 3c does not strictly follow the 
mean zonal wind contours, as meridional component is miss- 
ing here. 

There is no significant difference between two outmost 
views at the lOhPa winter pole and tropical lower strato- 
sphere, which implies the same amount of eastward propa- 
gating and westward propagating GWs. Simulation results 
from Sect. 4.3 confirm this inference. 

Very similar features can also be obtained with the 3-pt 
filter, though the amplitudes are smaller (Appendix B), con- 
firming that the result is robust. 


3.2 Geographical distributions 

Monthly mean geographical maps are plotted in Figs. 5 and 6 
at various levels to represent typical situations of January and 
July. Only the variance that is larger than 2er m j n is colored, 
where er m i n is the minimum detectable variance. Zonal mean 
wind (black solid contours) and wind speed at 700 hPa that 
exceeds 10 ms -1 (brown hatched areas) are overplotted on 
the variance maps. 

Like the zonal mean maps Figs. 3 and 4, GW variances 
show great enhancements over middle-high latitudes in the 
winter hemisphere where the polar night jet is the strongest. 
Clear orographic GW signals can be identified amongst the 
overall enhancement. In the northern hemisphere (NH), GW 
variances as well as the differences show great longitudinal 
dependence. In January 2005, a clear subtropical jet and a 
mid-latitude jet can be identified from the 700 hPa wind. The 
mid-latitude jet (which roughly collocates with the upper- 
level polar night jet) is apparently more powerful in gener- 
ating the MWs when it encounters big south-north oriented 
mountain ranges such as Rockies, the tip of the Appalachian 
mountains (Mts), Greenland, Iceland, and Kjolen Mts at the 
Scandinavia Peninsula. Since the polar night jet in the strato- 
sphere is also strong and favorably aligned with the tropo- 
sphere jet over Russia, enhancements can be seen over some 
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(d) West, 10 hPo x 10” (K 1 ) ( e ) East, 10 hPa 


(f) West - East, 10 hPa x io”(K’) 



(g) West, 80 hPo x 10 ” (K 1 ) (h) East, 80 hPa (i) West - Eost, 80 hPo x io”(k') 



Fig. 6. Same as Fig. 5, except for July 2005. 


Russian Mts as well. The subtropical jet, on the contrary, 
barely generates any large amplitude GWs. That is prob- 
ably because the subtropical jet center locates in the upper 
troposphere while the polar jet in the troposphere is closer 
to the ground, as can be found in Fig. 3. Also, note that 
wind above the subtropical jet switches direction in the mid- 
dle to upper stratosphere. One can imagine GWs generated 
in this region must suffer from critical level filtering or turn- 
ing point reflection (c.f., Dunkerton, 1984). In the Southern 
Hemisphere (SH), the two jets in the troposphere are more or 
less mixed together, and huge MWs are generated over the 
tip of the Andes, tip of the Antarctic Peninsula, as well as 
some islands (e.g., New Zealand). In the upper stratosphere, 
the Andes seem to have greater impact downstream as the 
emitted westward propagating GWs present across the globe 
at that latitude (top right plot in Fig. 6). 

Such orographic GWs have been seen through other satel- 
lite instruments, such as Aura MLS (Wu and Eckermann, 
2008) and AMSU-A (Wu et ah, 2006). The GW ampli- 
tudes in AIRS, however, are the strongest among all even 
after taking into account the averaging effect. This indicates 
that mountain waves (MWs) likely possess greater amount 
of energy toward the high-frequency part of the frequency 
spectrum (relative to the mean wind), and AIRS is particu- 
larly sensitive to them due to its high resolution and the steep 
viewing angle. In both hemispheres, orographic GW maxima 
are at high-latitudes (Greenland, Antarctic peninsula) in the 
lower stratosphere, and the peaks within 50-60 deg become 


more and more identified in the middle and upper strato- 
sphere (eastern Canada, south tip of Andes). This feature can 
be attributed to vertical wind structures, as the strong wind 
near the jet center Doppler-shifts the GWs toward longer ver- 
tical wavelength, and GWs near the core of the jet hence have 
better visibility to AIRS. The orographic GWs are closely as- 
sociated with the upper level polar night jet, and the location 
of the enhancements vary from year to year with the shift 
of the polar night jet. Comparing the lOhPa map with the 
2.5 hPa map, significant increase of GW activities is found 
in Alaska in January 2005 (Fig. 5) and West of Antarctic 
Peninsula in July, 2005 (Fig. 6). They are highly variable 
along the longitudes, and not necessarily related with the jets 
since they do not locate at the exist region of the jet stream, 
nor do they show preferred propagation directions. Inter- 
estingly, the Brunt-Vaisala frequency (N 2 ) also peaks at the 
same locations for these two particular months (not shown). 
In general, GW amplitudes are much smaller at low altitudes 
(80 hPa) than at higher levels. At high latitudes the GW ac- 
tivity decreases with altitude independent of longitude and 
not directly related to the wind velocities. This causes the 
high-latitude maxima around lOhPa in Figs. 3 and 4. We 
will discuss about this feature later on in Sect. 4.3. 

In the subtropics and tropics, large GW activities are found 
in the upper stratosphere over the deep convective regions. 
The deep convective regions are identified from the ice wa- 
ter content (IWC) from Aura MLS (Wu and Eckermann, 
2008). In particular, they are Western Pacific warm pool 
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Fig. 7. Multi-year averaged timeseries of zonal mean GW variance for different months at westmost view (first column), eastmost view 
(second column) and their difference (third column) at 2.5, 10, 40 and 80hPa from top to bottom. The UK Met Office zonal winds are 
overplotted in solid (dashed) lines representing westerlies (easterlies) with contour interval of 10 m s _ 1 . 


region, Amazon rainforest region, and Central Africa rainfor- 
est region for NH winters, and Southeastern US monsoon re- 
gion and India-South China monsoon region for SH winters. 
Comparing with Aura MLS, the magnitude of these convec- 
tively generated GWs seen from AIRS are smaller. This is 
consistent with the convective source spectrum suggested by 
Beres et al. (2004), where source momentum fluxes decrease 
rapidly at high-frequency part of the spectrum. The AIRS 
waves tend to propagate eastward relative to the mean wind, 
but the differences between the two outmost views have val- 
ues that are marginally above the minimum detectable ones, 
which means both eastward and westward propagating GWs 
exist in the upper stratosphere, and the westward propagat- 
ing ones are smaller than the eastward propagating ones with 
respect to the integrated magnitude. Convection is known 
to generate a broad spectrum of GWs that propagate in both 
directions at top of the deep convection (e.g., Alexander et 


al., 1997; Beres et al., 2004). Afterwards, the majority of 
the low-frequency westward propagating GWs are filtered 
by the easterlies in the tropics and subtropics, but some of 
the high-frequency waves may survive because they gener- 
ally have higher phase speeds. These convectively excited 
high-frequency GWs do not grow into a comparable strength 
with MWs until they reach the upper stratosphere. In the 
lower stratosphere, a imifonn belt of GW enhancements is 
found in the tropical region, with some highlights in the deep 
convective regions, which creates the tropical maximum in 
the lower stratosphere in the zonal mean maps. This belt has 
been seen by previous studies using various GW observa- 
tional instruments (e.g., radiosonde, GPS, IR limb sounders, 
etc., cf. Preusse et al., 2009). It apparently is not closely as- 
sociated with the convective sources. It is very interesting 
that AIRS can see it but AMSU-A, as another nadir viewing 
sounder, does not see it at all. The belt will be revisited later 
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on in Sect. 4.3, where we argue that the belt is caused by the 
propagation of inertia GWs there. 

In the high latitude summer hemisphere, GW signals are 
barely detectable. The 3-pt results without subtracting noise 
yield a quite uniform distribution of “background-like” vari- 
ance, which is independent of location and wind. In the NH 
summer (Fig. 6), it is easy to explain the low wave activity 
in the NH since the wind at 700 hPa is too weak to generate 
large orographic GWs (brown hatched regions barely exist in 
NH in Fig. 6). In the SH summer, we do see some enhance- 
ments over the Andes in the lower stratosphere (Fig. 5i), but 
GWs hardly propagate upward probably because of the criti- 
cal level filtering effect at the zero-wind line. 

Large amplitude MWs vary their location year to year (not 
shown). The 10ms _1 wind speed at 700 hPa roughly pro- 
vides a good threshold for generating noticable topographic 
GWs in the stratosphere. Only at places where westerly wind 
is consistently strong from low troposphere to upper strato- 
sphere can those GWs grow into significant strengths. Loca- 
tions of high occurrence of convectively generated GWs are 
also closely associated with the movement of tropical deep 
convective zones. 

3.3 Temporal variations 

A pronounced annual cycle with maximum in the winter 
hemisphere over mid-high latitudes is found throughout the 
region of interest in both the variance and difference time- 
series. We pick the levels of 2.5, 10, 40 and 80hPa for the 
timeseries in Fig. 7, together with the monthly zonal wind 
derived from UK Met Office dataset. Since the annual cy- 
cle is highly repeatable, multi-year AIRS variances are av- 
eraged together in Fig. 7 as well as the UKMO winds. The 
annual cycle found in AIRS are observed, for instance, by 
Aura-MLS (Wu and Eckermann, 2008), radiosonde (Wang 
and Geller, 2003), GPS (Tsuda et al., 2000), IR limb sounder 
(Preusse et al., 2009), etc. Eckermann (1995) attribute the an- 
nual cycle in both lidar and rocket sounding observations to 
seasonal variations in the density stratification of the atmo- 
sphere (i.e., N 2 ). Besides these major factors, source proper- 
ties can also contribute to the variations in that stronger west- 
erly winds during hemispheric winters tend to create larger 
MWs, which propagate further downstream under preferable 
conditions. 

The Southern Hemisphere differs from the Northern 
Hemisphere in terms of the location of the maximum vari- 
ance. In the Northern Hemisphere, wave variance nearly al- 
ways peaks at the time when jet is the strongest (December- 
Feburary), as do the differences. Interestingly, in the South- 
ern Hemisphere, the variance peak occurs ahead of the peak 
of the westerly jet below lOhPa, but the difference between 
the two outmost viewpoints follows the variation of the jet 
closely, both of which propagate northward in time (May 
through October). Below lOhPa, GW amplitudes are largest 
in July, but the mean wind maximum comes in one month 


later. A close look at the vertical structure of the zonal wind 
reveals that tropospheric westerlies in the SH high latitudes 
actually peak in July (not shown). Another plausible expla- 
nation is based on the fact that the coastline of Antarctica 
is approximately parallel to the latitude lines. Therefore, 
meridional winds are more important than zonal winds there 
in generating large amplitude waves. Actually, the merid- 
ional wind contours follow the change of GW variance more 
closely than zonal wind contours (not shown). The Andes 
seem to be vital at high levels (above lOhPa) at the intensi- 
fication stage of the jet speed, while Antarctic Peninsula is 
more important in lower levels and seems to be a major fac- 
tor causing the northward movement of the center of polar 
night jet. This phenomenon is also evident in Fig. 6. Note 
that the polar night jet in SH winter (e.g., Fig. 3) curves pole- 
ward (equatorward) below (above) lOhPa. Antarctic Penin- 
sula and the Andes perhaps play different roles at different 
levels in the formation of this shape. In the NH winter, there 
is no such northward curvature in either the jet center or the 
center of wave variance maximums. This result again sug- 
gests the importance of mountain GWD in shaping and reg- 
ulating the polar night jets. 

Strong interannual variability exists in NH high latitudes 
(not shown). The polar night jet splits during the winters of 
2004, 2006 and 2009 at the polar region, and GW variance 
reduces a significant amount simultaneously. These are the 
years of strong sudden stratospheric warming (SSW) events 
(Lee et al., 2009). Lee et al. (2011) observed enhanced plan- 
etary wave activities during the SSW events from Aura MLS 
observations, and they believe the planetary waves serve as 
an efficient filter to prevent upward propagation of GWs. 
More explanations can be found in Dunkerton and Butchart 
(1984). 

The annual cycle also dominates the subtropical strato- 
sphere, but completely opposite phases are found between 
the upper stratosphere and lower stratosphere. In the up- 
per stratosphere (first row of Fig. 7), GW variance as well 
as west-east difference peak during hemispheric summers. 
Zonal mean maps (Figs. 5 and 6) show that the convectively 
generated GWs grow significantly in the upper stratosphere, 
which collocate with the easterly jet centers as well. There- 
fore, this summer peak could be largely explained by the con- 
vective source properties. In the lower stratosphere, winter 
maximum dominates, and this should be again attributed to 
the variations of density stratification similar to that of the 
mid-high latitudes (Eckermann, 1995). 

The tropical lower stratosphere has an annual cycle with 
NH winter maximum, and this largely agrees with the fact 
that deep convections in the SH summers are in general 
stronger (Liu et al., 2007). QBO signals are relatively weak 
but still significant compared with the annual cycle, show- 
ing very interesting features that will be discussed in the 
next section. In the middle to upper stratosphere, a prevail- 
ing semi-annual cycle is found (not shown), which we be- 
lieve is a combined effect of background winds (semi-annual 
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Fig. 8. (a-c) Geographical map of GW variance at eastmost view at three different locations for January 2005 at 2.5 hPa. The GW variances 
within the white box of (a-c) as a function of FOV numbers are plotted in solid for each location at (d-f) January 2005, 2.5 hPa; (h-j) 
January 2008, 2.5 hPa; (k-m) January 2005, 80hPa; and (n-p) January 2008, 80 hPa. Negative (positive) FOV numbers correspond to west 
(east) views, and 0 corresponds to nadir view. The linear fits of each FOV curve are overplotted in dash-dotted lines. 


oscillations there) and the source (convection crossing the 
equator twice per year at equinoxes). 

4 Discussions 

Some new features from AIRS GW observations are worth 
further discussions. In this section, we will provide in- 
depth analyses on GW sources, interannual variations, and 
enhancements from inertia GWs. 

4.1 Orographic and convective GW sources 

Topographic GWs have been extensively studied during the 
last decades, and relatively mature parameterization schemes 
of the mountain GWD have already been employed widely in 
state-of-art GCMs. These schemes usually take into account 


factors such as mountain height, surface wind speed, atmo- 
spheric stability and the alignment of the mountain ridges 
(McFarlane, 1987; Bachmeister, 1993). The non-stationary 
convective GWs received less attention compared to the oro- 
graphic sources until recent years (Song and Chun, 2005; 
Beres et al., 2005). As to the parameterization of these 
GWs, some treat them as “moving mountains”, and some 
think they create much broader spectra (e.g., Song and Chun, 
2005). Both of the orographic and convective GWD param- 
eterizations have free parameters to represent the wave in- 
termittency, which is, by far, the poorest constrained param- 
eter. Other “tunable” parameters include horizontal wave- 
length, source momentum spectrum, etc. AIRS, as one 
powerful GW observation instrument, can certainly improve 
our global understanding on the “observational constraints” 
of the parameterization schemes. For example, from the 
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geographical distribution that has been discussed in Sect. 3.2, 
we can tell that MWs tend to have larger energy variance to- 
ward the high-frequency part of the energy spectrum, while 
it seems to be the other way around for the convective GWs. 
Besides, some interesting features are found in AIRS, which 
are quite different from what’s been observed before by other 
instruments. 

To look into more details, we pick two orographic GW 
cases and one convective GW case that are evident in January 
2005, which occur at the Rockies, the warmpool region, and 
Iceland. The monthly mean FOV-dependent GW variance 
within the white grid box is plotted at 2.5 hPa and 80 hPa for 
both January 2005 and 2008 in Fig. 8. At first glance, the lin- 
ear regression line confirms our idea that orographic GWs in 
general propagate westward while tropical convective GWs 
propagate eastward, both relatively to the local mean zonal 
wind. Ideally, if a monochromatic wave exists, we should 
see a smooth curve with a single peak or no peak depending 
on whether the upward propagation angle is smaller or larger 
than 48°. However, multiple peaks stand out in almost all 
curves in Fig.8d-p, indicating a spectrum of GWs are gen- 
erated down below. In the lower stratosphere (Fig. 8k-p), 
convectively excited GWs may have sporadic spikes (Fig. 8i, 
o) while orographic GWs can generate a rather broad spec- 
tram (Fig. 8n). The peaks of the variance in certain cases 
are close to the nadir (e.g., Fig. 8o, p), which means those 
Doppler-shifted GWs propagate almost vertically. By com- 
paring the 2.5 hPa FOV curves with the ones at 80 hPa, we 
may conclude that background wind plays a leading role in 
shaping the GW spectrum during the upward propagation. 
However, three preferred angles at FOV No. 25, No. 0 and 
No. 20 persist from lower to upper altitudes over the Rock- 
ies during January 2005, and the FOV No. 0 and No. 20 also 
have peaks of GW variance at 2.5 hPa during January 2008, 
even though no big wave event occurred over the Rockies 
in January 2008. Similar phenomena happen over Iceland 
(Fig. 8f,j), and the Andes (not shown). These preferred direc- 
tions should be associated with the detailed structures of the 
topography, and background wind is probably another candi- 
date that prefers to select certain angles. 

GWs seen by AIRS provide at least two new discoveries 
about the orographic and convective GWs. For orographic 
GWs, vertical wavelength is already large at the wave genera- 
tion level compared with convective GWs, and certain angles 
are favored for generating large amplitude GWs, which vary 
among different mountain ranges and vertical wind structures 
as well. The ratio between peaks and troughs as well as max- 
imum amplitude occurring FOV angles can be used to infer 
wave intermittency, which is an important wave parameter 
in GCMs that lacks observational constraints. It is not clear 
why certain angles are preferred. The reason why large am- 
plitude GWs can be seen in the nadir view is another inter- 
esting point and beyond the scope of this study. They remain 
as potential topics for future investigations. 


Convectively generated GWs, as suggested here in AIRS 
observations, have much smaller wave scales in terms of both 
vertical and horizontal wavelengths. Therefore, it is expected 
to have fair amount of momentum flux in the high-frequency 
end of the momentum spectrum, which also varies quite a lot 
based on sporadic spikes present in AIRS. The current as- 
sumptions of source spectrum (e.g., Song and Chun, 2005; 
Beres et al., 2004) for convective GWs is, at least, not accu- 
rate enough to represent all situations. Note that these curves 
are all derived from monthly means. Spikes in both moun- 
tain and convection cases may simply come from different 
events, i.e., wave intermittency is a highly variable parame- 
ter. More exploration of the details need to be carried out in 
the future before advice can be given to modelers. 

4.2 Interannual variations at the equatorial region 

The quasi-biennual oscillation (QBO) wind acts as an ef- 
ficient GW filter in the tropical lower-middle stratosphere, 
and the interations between GWs and background winds 
are, in turn, an essential ingredient in the QBO. Different 
GWs contribute differently to driving the QBO. Kawatani 
et al. (2010a, b) were able to simulate a realistic QBO 
with a high-resolution climate model without parameterizing 
GWD. Their GCM can resolve part of the inertia GW spec- 
trum (limited by the model’s vertical grid resolution), and 
the simulated QBO has a shorter period than that of the real 
atmosphere, but with realistic strength and structure. Their 
work suggests the importance of inertia GWs in generating 
the QBO. Some other model studies suggest that the momen- 
tum deposition associated with the filtered GWs accounts for 
at least half of the total momentum flux in driving the QBO 
(Giorgetta et al., 2002; Kawatani et al., 2010a). QBO signals 
have also been seen in various GW observations, such as ra- 
diosonde (Vincent and Alexander, 2000; Wang and Geller, 
2003), GPS (Torre et al., 2006), IR limb sounders (Krebs- 
bach and Preusse, 2007) and Aura MLS (Wu and Eckermann, 
2008). These instruments are ideal for observing low to mid- 
dle frequency GWs. 

AIRS GWs contain a strong annual cycle in the tropi- 
cal lower stratosphere (Fig. 7). On top of the annual cy- 
cle, a clear QBO signal stands out. After removing the an- 
nual cycle, the linear trend, and sub-seasonal oscillations, we 
present the timeseries of the GW variance/difference anoma- 
lies at the equator in Fig. 9. Pronounced GW variance en- 
hancements occur along with the descent of the QBO east- 
erly phase (defined as when the monthly mean zonal wind 
is easterly, likewise for QBO westerly phase), and a reduc- 
tion of GW activity happens right before and lasts through 
the whole period of the descending of the westerly phase 
(Fig. 9a). Year 2003-2004 is somewhat different from other 
years, when the enhancement only occurs during the down- 
ward propagation period of QBO easterly shear, and a strong 
negative anomaly occupies the rest of the phase. The GW 
differences also show preferred zonal propagation directions 
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Fig. 9. Monthly mean timeseries (annual cycle and linear trend re- 
moved) of the zonally averaged GW variance taken at the westmost 
view (a) and the variance differences between the westmost and the 
eastmost view (b). The vertical cross-section is taken at the equa- 
tor from 80 hPa to 20 hPa. 3-month running smooth window is ap- 
plied to the data to remove the sub-seasonal cycle. UKMO monthly 
mean zonal winds are overplotted in lines with contour intervals 
of 5ms -1 and solid (dotted) lines indicating eastward (westward) 
winds. Variance within ±0.085 x 10“ 1 K" and difference within 
±0.033 x 10 -3 K 3 are uncolored. 


during different QBO phases. Although it is more ambiguous 
in Fig. 9b, in general the GWs that can survive the wind per- 
fer to propagate eastward (westward) in the presence of QBO 
easterly (westerly) shears. Again, year 2003-2004 is an ex- 
ception. The other thing to notice is that the QBO can only 
account for about 1/10 of the variance that can be explained 
by the annual cycle. 

The close correlation shown in Fig. 9 indicates the in- 
teractions between QBO winds and some of the AIRS ob- 
served GWs. Moreover, it suggests that these GWs may 
play a more important role for the descent of QBO west- 
erly phase than that of the easterly phase based on the ev- 
idence that the suppression (enhancement) of wave activi- 
ties occurs at the downward propagation of QBO westerly 
(easterly) shear, i.e., more GWs are removed and hence de- 
posit their momentum fluxes at the downward propagation 
stage of QBO westerly shear. Deep convection is the most 
dominant source at the tropics that is responsible for GWs 
we see in AIRS, as discussed in the last section. Since the 
upper-troposphere wind is in general easterly, the spectrum 
of deep-convection generated GWs at the tropopause level 
inherently has a skewness toward positive phase speed, i.e., 
eastward propagating GWs dominate. Hence, QBO west- 
erlies tend to filter out more GWs than QBO easterlies, re- 
sulting in the observed less GW activity during the QBO 
westerly phase. This finding is consistent with what has 
been observed by Aura MLS (Wu and Eckermann, 2006). 
Some of the previous studies using radiosonde and GPS also 
found suppression of GWs in the QBO westerly phase (e.g., 


Sato and Dunkerton, 1997; Wang and Geller, 2003; Torre et 
al., 2006), but they suggest that this suppression occur right 
along the descending of the zero-wind line of the QBO west- 
erly phase. This is mainly caused by Kelvin waves that are 
unavoidably blended in the measurements from radiosonde 
and GPS. MLS scans along-track (almost meridionally over 
most of the latitudes), and Kelvin waves are truncated by 
a technique with short horizontal wavelength cut-off. 7-pt 
running smooth window we applied on AIRS measurements 
produces a cut-off wavelength at ~ 105 km, which can also 
effectively filter out Kelvin waves. Another possible reason 
lies in the fact that radiosonde and GPS observes more low- 
frequency GWs than AIRS does, which are believed to have 
larger impact on the QBO formation and downward propa- 
gation than high-frequency GWs (Kawatani et al., 2010a). 
Vincent and Alexander (2000) found increased GW activ- 
ities during QBO easterly phase through radiosonde data 
over a tropical island, which is a bit different from other ra- 
diosonde observations but is consistent with what we observe 
in AIRS. The reason for the strong negative anomaly during 
year 2003-2004 remains unclear. It might be associated with 
other interannual variabilities such as El Nino Sothern Oscil- 
lation (ENSO). Due to the limit of AIRS data record length, 
this conclusion is not definitive, and remain a nice topic for 
future investigation. 

However, the QBO only explains about 10% of the total 
variability of AIRS GWs, which means AIRS GWs are only 
slightly modified by the QBO rather than playing a domi- 
nant role on the formation/propagation of QBO phases. This 
is expected as AIRS is mostly sensitive to higher-frequency 
GWs, which usually propagate into the upper stratosphere 
and above (Fritts and Alexander, 2003). In fact, the max- 
imum wind speeds of the easterly and westerly associated 
with the QBO are —35 ms -1 and 20ms -1 . At the equator, 
this converts to a vertical wavelength range of 5.6 to 10 km 
assuming N 2 — 5x 10 4 s -2 . These GWs are too short to be 
seen by AIRS. Nevertheless, AIRS GWs still contain QBO 
signals, indicating that AIRS should have some capability of 
observing relatively short vertical wavelength waves, which 
will be explained in the next section. 

Theoretical work predicts both westward and eastward 
propagating GWs should both contribute to the QBO for- 
mation (Baldwin et al., 2001). This is further confirmed by 
simulation results from high-resolution GCMs (e.g., Sato et 
al., 1999; Kawatani et al., 2010a, b). However, it is difficult 
to separate eastward propagating GWs from massive Kelvin 
wave signals in radiosonde and GPS observations. The in- 
ferred preferred zonal wave propagation direction from AIRS 
at the equator show both eastward and westward propagation 
signals with about equal probability to occur (Fig. 9b). The 
uniqueness of AIRS GW response to QBO winds indicates 
the power of satellite instruments in obtaining the GW infor- 
mation. 
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Fig. 10. Zonal mean FOV curves at 80hPa for 3-pt (red) and 7-pt (black) results at January 2005 at 60° N (a), 0° (b) and 60° S (c). The 
horizontal axes are the number of scanning points along a FOV scan, where — (+)45 corresponds to westmost(eastmost) view, while 0 is at 
nadir. 


4.3 Wave enhancement in the equatorial lower 

stratosphere and near the winter pole at 10 hPa 

MWs are excessively strong in the AIRS observations, which 
overwhelm other GW signals if we look globally. In the 
lower stratosphere, however, GWs are stronger in the equato- 
rial region. These GWs are not closely associated with deep 
convection, as the variance forms a uniform belt (Figs. 5g, 
h, 6g,h, and bottom row of Fig. 7). By revisiting Fig. 3, 
we can find that the amplitude of these GWs maximizes 
at the tropopause (~ 100 hPa), decreases dramatically right 
above, and increases again above 80 hPa. These features are 
in good agreement with previous studies using radiosonde 
(e.g., Wang and Geller, 2003), GPS (e.g., Ratnam et al., 
2004), Aura MLS (Wu and Eckermann, 2008), CRISTA and 
SABER (Preusse et al., 2006), and high-resolution GCM 
simulations (e.g., Sato et al., 1999), but are not shown in 
UARS MLS (Wu and Waters, 1 996) nor Aura AMSU-A (Wu, 
2004). Meanwhile, there is also a peak at the winter pole at 
lOhPa, as shown in Figs. 3 and 4. This is the first observa- 
tion of such a lOhPa winter-pole-enhancement. Since this 
feature is not uniformly distributed among the high-latitude 
hemisphere winters (Figs. 5 and 6), and the locations vary 
from year to year (not shown), it is believed not purely caused 
by underestimation of the instrument noise, but rather a real 
atmospheric phenomenon. 

A careful evaluation of FOV-dependent variance profile 
(“FOV curve” hereafter) suggests that the 80 and 10 hPa vari- 
ance enhancements are something other than high-frequency 
internal GWs and instrumental noise. The zonal mean FOV 
curves derived at 60° N, 0° and 60° S for January, 2005 at 
80 hPa are plotted in Fig. 10. At 60° N, as we would expect to 
see from topographically generated internal GWs in North- 
ern Hemisphere winter, the variance is the largest at the east- 
most view. At 60° S, the FOV curves are basically flat and 
small, close to the estimated noise level. At the equator, the 
FOV curve bends downward, with maximum value occurring 
at the nadir view. This is also true at the two flanks of the po- 
lar night jet at middle to upper stratosphere (not shown). 

Turbulence can cause “bell” shape FOV curves seen in 
Fig. 10b due to the fact that the power index for 2-D tur- 
bulence is more negative at the off-nadir view than the nadir 


view and hence the area integral is smaller for the off-nadir 
view (Gruninger et al., 1998). Since shear instabilities eas- 
ily occur at the two flanks of the polar night jet, turbulence 
can be particularly strong there. However, we cannot explain 
why turbulence is particularly strong in the equatorial lower 
stratosphere and the winter pole at lOhPa with this theory. 
Besides, according to Aumann and Miller (1994), the out- 
most view has a FOV swath area 4 times larger than that of 
the nadir view, which would result in a very deeply down- 
ward curved FOV-dependent variance profile. Even with 
noise being already subtracted, the FOV curve in Fig. 10b 
has a difference between outmost and nadir views of only 
7 %, which is way too small to be purely interpreted as the 
turbulence. 

Another explanation of the downward curvature of FOV 
curves in Fig. 10b invokes short vertical wavelength, long 
horizontal wavelength, low-frequency inertia GWs. These 
GWs are probably the cause of the QBO signal presented 
in the last section. Alexander et al. (2002) proposed the 
idea that slow GWs have larger probabilities to be observed 
than fast waves. Since the vertical group velocity is approxi- 
mately proportional to square of the wave intrinsic frequency, 
AIRS has a larger chance to observe inertia GWs in the lower 
stratosphere, and the overall enhancement of wave activity 
observed in the tropics is largely due to the reason that these 
inertia GWs propagate slantwise and spread out from the 
convective source consequently. Since the GW enhancement 
in the equatorial lower stratosphere observed by AIRS is sim- 
ilar to GWs revealed by radiosonde and GPS as well as sim- 
ulation results in high-resolution GCMs (e.g., Fig. 8c in Sato 
et al., 1999), we have more confidence to claim that AIRS 
observes inertia GW signals. 

Here we conduct a simulation to show how AIRS detects 
inertia GWs of short vertical wavelength in the lower strato- 
sphere. Assuming a constant wave amplitude, only waves 
with vertical wavelength longer than 12 km are visible in 
the nadir viewing center tracks. The WF is practically a 
delta function with respect to long vertical wavelength GWs 
and therefore yield GW variance at the level where the WF 
peaks. For waves with vertical wavelength less than 12 km, 
the convolution of waves of a constant amplitude with AIRS 
WF will yield a very small variance. If, however, the wave 
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Fig. 11. (a) Weighting function over nadir at 80 hPa as a function of 
height and along-track distance for the numerical simulation. The 
contours are 0.5, 0.7, 0.9 and 0.99 of the peak response. The out- 
most off-nadir weighting function is 4 times broader than that of the 
nadir, (b) One example of the original input wave convolved with 
the nadir weighting function (lines) and resultant wave with 0.15 jt 
phase shift (shades). The wave has horizontal wavelength (A*-) of 
100 km and vertical wavelength (X z ) of 6 km. Negative values are 
contoured with dash-dot lines (light colors), and positive values are 
denoted with solid lines (dark colors). The contour interval is 0.2 of 
the maximum and minimum amplitudes, respectively. 


amplitude varies with height, this can result in an effective 
narrowing of the weighting function, i.e., the GWs at the 
edge of the WF cannot be completely cancelled, and will 
yield some response in AIRS variance. Figure 11a gives an 
ideal WF at nadir view at 80hPa to represent mean WFs of 
multiple AIRS channels we select which all peak at 80 hPa 
(Appendix A). The WF at the outmost off-nadir view cov- 
ers an area 4 times broader than that of the nadir WF (see 
Appendix C for details). Figure lib shows the GW that is a 
result of the original imported GW convolved with the WF. 
The original GW has an amplitude varying with N 2 , and the 
maximum amplitude A max is 5 K (see Appendix C for de- 
tails). 

The final goal of this experiment is to try to find the 
best wave parameters that give the convolved variance at the 
nadir and the variance difference between the nadir and out- 
most views consistent with the observations (black line in 
Fig. 10b). In order to compute the variance, we make a hor- 
izontal phase shift of nrc forward of the convolved wave, as 
shown in shades of Fig. lib. This procedure is to account 
for the 7-pt smoothing window we applied to AIRS data. 
The difference between the two convolved waves are then 
calculated and compared with the observation (Fig. 10b) at 
the nadir view. With a fixed X x of 700 km, 7-pt smoothing 
window length (~ 105 km) corresponds to 0.157T phase shift. 
We can therefore plot out the expected AIRS variance as a 
function of vertical wavelength X z in Fig. 12a. The same 
procedure can be applied to the outmost off-nadir footprint 
to obtain a reading of wave variance there, and the differ- 


Fig. 12. Variance differences between unshifted wave and the wave 
shifted by 0. 15 jt for X x = 700 km calculation (a), and the variance 
difference between nadir FOV and the outmost FOV (b) as func- 
tions of X z . The AIRS observed values (7-pt smoothing window 
applied) are denoted by dash lines. See text for details. 


ence between the nadir and outmost views can also be com- 
puted as shown in Fig. 12b. The two curves meet with the 
observed values both at X z — 6 km. We tried other X x val- 
ues and it turns out X x — 700 km, X z — 6 km and A max = 5 K 
gives the best agreements simultaneously. Typical horizon- 
tal and vertical wavelengths for equatorial inertia GWs are 
around 1000 km and 5 km, respectively (Wang and Geller, 
2003), and the peak amplitude is around 5K (Sato et al., 
1999), which is in good agreement with our simulation re- 
sult. The inertia GW variance can also explain why there is 
no obvious propagation direction in the difference between 
eastmost and westmost views, as the radiance response in- 
volves no propagation direction information whatsoever in 
this case. The layer near 80 hPa (Fig. 1 la) is less affected by 
the QBO wind than the layers above and hence we still do 
not see a dominant QBO signal even though these waves are 
believed to be dominantly inertia GWs. The simulation with 
a realistic AIRS WF supports the idea that the enhanced GW 
variance at lOhPa winter pole is also low-frequency inertia- 
gravity waves. See Appendix C for details about this experi- 
ment. 

There are quite a few explanations of why the equatorial 
maximum weaken rapidly with height. One possible reason 
is that these inertia GWs with low vertical group velocities 
simply dissipate. If they speed-up by refraction, the visibil- 
ity may lessen. Another possibility is that these waves meet 
the “tropical confinement”. Once they propagate off-equator 
due to the slantwise propagation, they tend to easily meet the 
lower limit of the frequency, / (inertial frequency), as / in- 
creases with latitude (Sato et al., 1999). Hence the tropical 
belt of GW enhancement disappears at high altitudes (20 hPa 
and above), and only those high-frequency internal GWs that 
are closely related with convective sources can survive. Choi 
et al. (2011) applied 3-D AIRS WFs to their parameterized 
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convective GWs, and found out that merely any GWs can 
be seen in the equatorial lower stratosphere since the param- 
eterized GWs are too small in terms of the vertical wave- 
length and horizontal wavelength as well. Propagating up- 
ward, these small GWs are Doppler-shifted toward longer 
vertical wavelength, and become detectable by AIRS, which 
agree very well with what we see at middle and upper strato- 
sphere. 

A peak at 10 hPa at the winter pole appears in the N 2 field 
in Fig. Cl of the Appendix C. Therefore, as in the tropical 
lower stratosphere, inertia GWs should grow again at those 
levels, but with a short penetration. Convolving with the 
weighting function, and the GWs result in another enhance- 
ment in AIRS variance. Due to the sensitivity limitations, ra- 
diosonde and GPS sensors cannot observe those GWs at this 
level. UARS MLS and Aura AMSU-A cannot see these fea- 
tures probably because of the coarse horizontal resolutions 
(Wu and Eckermann, 2008). Aura MLS also observes sim- 
ilar equatorial enhancement at lower levels, but not at the 
lOhPa winter pole, which remains to be investigated in the 
future (Wu and Eckermann, 2008). 

After all, the AIRS GW variance is believed to be com- 
posed of large amplitude internal GWs that are closely re- 
lated with the local sources (e.g., orography, deep convec- 
tions), and small amplitude “background-like” inertia GWs 
that are determined by the shape of the WFs and the shape of 
the original amplitude spectra. It is shown that AIRS can ob- 
serve the inertia GWs even though their vertical wavelengths 
are smaller than AIRS WF thickness. This interpretation 
might be applicable to explaining signals discovered in other 
instrument measurements. The two groups of GWs are po- 
tentially separable by regressing the variance on N 2 . Turbu- 
lence is another possibility, but the strength and mechanism 
remain as a question for future investigation. 

5 Conclusions 

In this study, climatology of the GW variance and zonal pre- 
ferred propagation direction measured by Aqua-AIRS are 
documented. In summary, evidence is found for short hor- 
izontal wavelength internal GWs (A./, ~ 100km, X z > 12 km) 
forming dedicated geographical maxima that are closely re- 
lated to localized sources such as topography and deep con- 
vection. A generally smaller background of longer horizon- 
tal wavelength GWs forms band-like structures for instance 
in the lower stratosphere tropics. It is remarkable that AIRS 
can detect these inertia GWs even though their vertical wave- 
length is smaller than the width of the AIRS weighting func- 
tion, if the wave amplitude changes notably over the range of 
detection. 

GW induced temperature variance is first succesfully ex- 
tracted from the brightness temperature variance observed 
by AIRS by carefully removing instrumental noises at var- 
ious pressure levels. A total of 50 AIRS channels is used 


with WFs peaking at various levels covering the entire strato- 
sphere. Because of the nadir viewing geometry and small 
footprint size, AIRS radiance variance is most sensitive to 
high-frequency internal GWs with comparable vertical and 
horizontal wavelengths with wave fronts in parallel to the 
scanning angle. Large smearing occurs when the wave front 
is across the AIRS scan, leading to possibility of inferring 
the preferred GW propagation direction using variance dif- 
ference between two outmost views. 

Comparisons have been carried out throughout the paper 
between AIRS GW variance and other measurements includ- 
ing satellite and non-satellite instruments (e.g., Aura MLS, 
Aqua AMSU-A, GPS, radiosonde, GCM simulations, etc.). 
The GW variance seen by AIRS in general agrees with pre- 
vious studies in terms of the distributions and variations of 
maxima. The AIRS GW peaks are generally found to be 
above high meridionally-oriented mountain ranges, as well 
as near tropical deep convection regions. Different from oth- 
ers, the GWs observed in AIRS are highly localized, since 
they are dominated by high-frequency components. This 
provides better information on GW sources from direct re- 
trievals of GWs, as shown in Alexander and Barnet (2007). 
The locations of GW peaks vary significantly from year to 
year. The 700 hPa wind is used to evaluate the GW gen- 
eration, showing that stratospheric GWs tend to grow into 
great amplitudes at the places with wind speed greater than 
10 ms -1 collocated with the upper-level jet maxima. In the 
tropics, the deep convection controls the GW variance peaks. 

The inferred preferred zonal propagation directions of 
GWs seen in AIRS are opposite to the mean zonal wind di- 
rections. Topographic GWs in the winter high-latitudes tend 
to propagate westward relative to the mean flow, while con- 
vective ly excited GWs in the tropics and subtropics tend to 
propagate eastward relative to the mean winds. The preferred 
wave propagation becomes clearer with increasing altitude, 
suggesting effective wave filtering and increasing GWD on 
the mean flow. The GW variance contours do not exactly fol- 
low the wind contours, suggesting less filtering of the high- 
frequency GWs, which propagate through and become more 
important in the mesosphere and above. 

An analogous approach of determining wave propagation 
direction has also been carried out on Aura MLS obser- 
vations (Wu and Eckermann, 2008), where wave preferred 
meridional propagation direction is inferred in a similar way 
considering the fact that MLS scans along the orbit. One can 
therefore obtain a more complete picture of the GW horizon- 
tal propagation direction from the combined AIRS and MLS 
results. Recent studies show that convective GWD param- 
eterizations are sensitive to the wave propagation direction, 
which is set to be a free parameter in the model (Song and 
Chun, 2005, 2008). Observations of the propagation direc- 
tion by AIRS may help to constrain this parameter. 

The new findings with AIRS GW variance differ from 
other GW observations in three aspects: 
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Table Al. Comparison of sensitivity to GW parameters from various satellite measurement techniques. Spectral coverage and sensitivity 
peaks depend on the techniques that are used in the corresponding reference. 


Instrument Name 

Altitude Range 

Spectral coverage 

Peak Sensitivity 

Reference 

AIRS 

surface —40 km 

X z >12km, X/, >13km 

X z = A/,~100km 

this paper 

GPS 

10-30 km 

X/j >100 km 

^■h /^z ^ 1 0 

Alexander et al. (2008) 

Aura-MLS 

15-50 km 

X z >5 km, A/, >50 km 

X z = 10 km, X/j = 100 km 

Wu and Eckennann (2008) 

CRISTA 

15-60 km 

X z >5 km, Xfr > 100 km 

x h /x z = m 

Preusse et al. (2009) 


Firstly, AIRS 90-footprint FOV-depedent GW variance 
above mountain ranges and deep convection occurs with pre- 
ferred angles. Because of the AIRS steep viewing angle, this 
finding indicates that really high-frequency GWs with large 
amplitudes are present, and mainly excited from orographic 
and convective sources. It requires further investigation on 
why certain angles are preferred, which is potentially a pa- 
rameter missing or poorly constrained (e.g., wave intermit- 
tency) in all GWD parameterization schemes. 

Secondly, the QBO signals observed in AIRS are different 
from other measurements that mainly observe inertia GWs 
and Kelvin waves. With relatively fast vertical group ve- 
locity, these high-frequency internal GWs seen in AIRS are 
more difficult to be filtered out by QBO winds, and hence 
they are less important in the formation and descending of 
QBO phases. AIRS GWs that are affected by QBO might 
be low-frequency inertia GWs that directly cause the lower 
stratosphere maxima of wave intensity. The QBO has not 
been well simulated so far in GCMs. On one hand, the coarse 
vertical resolution results in deficiency in resolving some 
inertia GWs. On the other hand, the poorly observation- 
ally constrained GWD parameterization schemes cannot pre- 
cisely represent the characteristics of high-frequency GWs. 
AIRS not only provides abundant information on the latter 
on to improve our understanding of high-frequency GWs, but 
also show potential to observe low-frequency inertia GWs. 

Lastly, the most important and encouraging finding from 
this work the AIRS capability of detecting short vertical 
wavelength inertia GWs, by assuming inertia GWs’ ampli- 
tudes varying only with background stabilities. AIRS ap- 
pears to be highly selective on high-frequency GWs, and 
therefore we literally can separate the GW variance between 
the inertia-gravity and high-frequency components. More- 
over, this new capability will help us to further understand 
GW characteristics as observed by other sensors. 

This paper provides the first climatological report of GW 
characteristics from nadir viewing sounders. Since AIRS is 
sensitive to high-frequency GWs, the results from this study 
also provide the first global survey of the characteristics and 
properties of high-frequency GWs. Some large-amplitude 
GW events as well as FOV-dependent variance can be stud- 
ied in further detail with models. Improving the GWD pa- 
rameterizations with the AIRS GW observations is another 
important task to explore in future study. 


Appendix A 

Channel numbers and the estimated noise level 

The channel numbers we used in this paper are listed be- 
low in Table A2. Together we also give the estimated noise 
from “3-pt running average window” method and the pro- 
vided instrument NEdT on 30 August 2002. The minimum 
detectable GW variances for zonal means and monthly maps 
are listed in the last two columns. Bold numbers in the 2nd 
column are channels that locate at the wings of the radiance 
spectrum and have no paired channels. These channels are 
particularly sensitive to small drift in the wavelength while 
others are stable. 


Appendix B 
3-pt analyses result 

As mentioned in the main text, 7-pt is the best filter size 
among the ones we have tested (3-pt, 7-pt and 15-pt). To 
better illustrate this point, 3-pt analysis results are shown in 
Fig. B1 and Fig. B2 (January and July of 2005, respectively). 
The overall structure is nearly identical to 7-pt results, but 
the amplitude is in general smaller than that of 7-pt analysis. 
We hence can conclude that the signals are robust, and GWs 
captured by 7-pt filter have larger power on the spectrum than 
those from 3-pt filter. 


Appendix C 

Some details about the simulation 


In the simulation mentioned in Sect. 4.3, we assume the im- 
ported wave amplitude varies with N. The reason is as fol- 
lows. Eckermann (1995) derived temperature variance as: 


-j- N z r 

T ' 2 =~ y 1 
8 L 


(f/cvY 


^M /2 + v' 2 ^ 


(Cl) 


where g, u', and v' are the gravitional acceleration, zonal 
wind perturbation, and meridional wind perturbation caused 

by GWs, respectively. T' 2 — T' 2 / 7' ( 2 , where T' is the 
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Table A2. A list of channel numbers, noise, NEdT, and minimum detectable GW variance at each pressure level. Bold numbers in the 2nd 
column are channels that are unpaired. See context for details. 


Min. detectable GW 
var. (x10“ 3 K 2 ) 


Pressure 

(hPa) 

Channel numbers 

Noise 

(K 2 ) 

NEdT 

(K 2 ) 

Zonal 

mean 

Map 

2 

74 

0.149 

0.165 

3.78 

26.64 

2.5 

75 

0.147 

0.166 

3.72 

26.22 

3 

76 

0.143 

0.161 

3.63 

25.55 

4 

77 

0.145 

0.160 

3.66 

25.80 

7 

78 

0.153 

0.162 

3.88 

27.34 

10 

79 

0.182 

0.172 

4.62 

32.53 

20 

81,82 

0.084 

0.078 

2.14 

15.05 

30 

102, 108, 114, 120, 125, 126 

0.039 

0.029 

0.98 

6.88 

40 

64, 88, 90, 94, 100, 106, 118 

0.033 

0.028 

0.83 

5.86 

60 

66, 68, 70, 86, 87,91,93,97, 130 

0.026 

0.018 

0.66 

4.68 

80 

92, 98, 104, 105, 110, 111, 116, 
117, 122, 123, 128, 129, 134, 140 

0.020 

0.011 

0.50 

3.54 

100 

132, 133, 138, 139, 149, 152 

0.026 

0.014 

0.67 

4.73 


(o) West (b) Eost 
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Fig. Bl. Same with Fig. 3, but 3-pt filter is applied instead of 7-pt filter. 


absolute temperature perturbation, while 7o is the back- 
ground temperature. The overbars denote the column aver- 
age. 

Think of it physically. Larger N 2 corresponds to more sta- 
ble atmosphere, which means the potential temperature in- 
creases more rapidly with height. Therefore, same oscilla- 
tions in height would cause larger temperature perturbations. 
By assuming that waves neither break nor meet saturations, 
the vertical group velocity varies inversely with N 2 , so that 
the wave activity coming up from below accumulates in the 
lowermost tropical stratosphere even if the waves are initially 
conservative there. 


The zonal mean values of N 2 at January, 2005 are shown 
in Fig. Cl. The input 2-D wave amplitude A is assumed to 
be proportional to N at the equator, which is approximated 
as: 


A(x,z) 


ry ry ^ ^60 ^ 

A 0 cos(^ + ^)e A , if z > z 60 

ry ry ~U~^60 I 

A 0 cos(^ + t f^)e , if z < Z60 

A z A x 


(C2) 


where Z60 ~ 21.5 km roughly corresponds to the geometry 
height at 60 hPa. x is the along-track distance, and z is the 
height. 
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Fig. B2. Same with Fig. Bl, except for 2005.07. 
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Fig. Cl. Zonal mean of the Brunt-Vaisala frequency (N~) as a 
function of latitude and height for January 2005, derived from ERA 
Interim dataset. 


An example of this weighting function has been included in 
Fig. 11a. At the outmost off-nadir view, the width of this 
weighting function on along-track direction is broadened by 
a factor of 5. 
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